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Summary. — The Grad-Shafranov approach to the problem of the structure of the 
black hole magnetosphere is discussed. For the double transonic flow, the number of 
boundary conditions in the pair creation region is shown to be sufficient to determine 
not only the longitudinal electric current, but also the angular velocity of a flow as 
a solution to a problem. As a result, the energy loss is determined by the physical 
parameters at the particle creation region rather than the "boundary conditions" at 
the event horizon. 

PACS 04.20 - Classical general relativity. 



1. Introduction 

The Blandford-Znajek [1] process of electromagnetic energy extraction from a rotating 
black hole is considered as the most appreciable source of activity of the central engine in 
Active Galactic Nuclei (see e.g. [2]). Indeed, the energy loss of a black hole with a mass 
m ~ lO^mQ and radius rjj — 2Gm/c^ ^ 10^^ cm embedded into an external poloidal 
magnetic field B ~ 10^ G and rotating with the angular velocity JIh [3] 




is large enough to explain the energy of jets and their radiation. Here < flu is 
the angular velocity of plasma, -BEdd ^ 10^ G is the standard estimate of the poloidal 
magnetic field near a supermassive black hole, and the factor fc '--^ 1 depends on the 
geometry of the magnetic field. The same process may play the leading role in some 
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galactic sources containing black holes of solar mass and even in the sources of gamma- 
ray bursts. 

In what follows wc use the standard expressions for the metric in Boyer-Lindquist 
coordinates (for c = G = 1) 

(2) ds^ = -a^dt^ + Qikidx' + I3'dt){dx^ + p''dt), 
where 

(3) a=|VA, P^ = f3' = 0, /3v = -a; = -^, 

2 2 
Qrr = gee = P , 9ipif = ^ ■ 

Here a is the lapse function vanishing at the event horizon 



(4) rn = TO + -\/to^ — a^, 
w is the Lense-Thirring angular velocity, and 

(5) A = + - 2TOr, = + cos^ 9, 

= (r2 + - A sin^ 9, w=- sin 9. 

P 

As a result, a/m = 2riH''H, where flu = ^'^(''h) is the angular velocity of a rotating 
black hole. In what follows, all three-dimensional vectors correspond to the 3 + 1 split 
language [3]. 

The Blandford Znajck mechanism of energy extraction is analogous in many respects 
to the one working in the pulsar magnetosphere [4, 5]. Indeed, in the presence of the 
longitudinal electric current / producing the toroidal magnetic field 

(6) ^^ = -7r 

aw 

and the angular rotation of plasma fip resulting from the electric field 

(7) E = -i^V*, 

lira 

the flux of the Poynting vector at large distances 

(8) Wem = ^ / [E X B]dS 

just corresponds to the value T4^bz (!)• Here ^{r,9) is the magnetic flux function deter- 
mining the poloidal magnetic field 
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It is necessary to stress the key point - we propose that in the magnetosphere there is 
enough plasma to screen the longitudinal electric field. Only in this case can the electric 
field E be presented in a form (7). 

As a result, the energy loss for a magnetically dominated flow can be rewritten in the 
form 



(10) Wbz= / 

Jo 



fiF(4')/(4') 



27r 

In other words, we have Wbz — IV^ax, where 

(11) Knax « ^^F*max/C 

is the potential drop between the central (^' = 0) and the marginal (^' = 5'max) magnetic 
surfaces passing through the horizon. The nonrelativistic relations (8) and (11) (in which 
we restore the dimension, c ^ I) can be used for the black hole magnetosphere as well 
because at distances r ^ th the General Relativity effects are insignificant. 

On the other hand, the effects of General Relativity were actually used in expression 
(1). The point is that the natural regularity condition at the horizon (the absence of 
infinite electromagnetic fields in the reference frame comoving with freely falling observer) 
Eg — > {Eg + B^)/a < 00, i.e. Eg + can be rewritten in the form of a "boundary 

condition" [6] 

(12) 4./w = i»H-»KW] ^/;Xt, .i.^(f 

It is the proportionality / oc (Oh — Op) that is responsible for the appropriate factor in 
(1)- 

Moreover, to forward the analogy, within the membrane paradigm it is possible to 
introduce a finite "surface charge density" (Th and "surface currents" Jh [3] 

(13) .H = 

(14) AnJu X n = B^e^ = aB^e^. 

As a result, the energy loss can be formally rewritten in the form 

(15) Wtot = j [EhJh - Phb^ (<JhEh + Jh X Bh)] dS, 

where /3h = P{i"h) and En = aEg. This expression is similar to the energy loss of radio 
pulsars 

(16) VFtot = "K. 
Here 

(17) K = y^[r X [Js X B,]]dS 
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is the braking torque resulting from the Ampere action of the surface currents Jg that 
close the electric currents flowing in the pulsar magnetosphere. Hence, one could imagine 
that it is the surface currents flowing along the horizon that result in the black hole 
braking. This physically clear picture was used later in many papers devoted to the 
central engine in Active Galactic Nuclei. 

Nevertheless, this approach met with some problems. First of all, in the last several 
years some critical papers appeared concerning the efficiency of the BZ process. The 
point is that, as one can see from (1), the necessary energy loss 10^^ erg/s can only be 
achieved for extreme values of the parameters. Actually, all parts of this equation were 
criticised. First of all, the Eddington magnetic field [2] 

(whose energy density is equal to that of the accreting matter giving Eddington lumi- 
nosity) is actually the upper limit of the magnetic field which can be generated by the 
accreting plasma in the vicinity of the black hole horizon. Up to now, the possibility of 
the generation of such a large regular magnetic field has not been confirmed by direct 
calculations [7]. Further, an extreme rotation a/m w 1 as well as very large masses 
m ~ 10^111,0 of a black hole arc necessary to reach an energy Wtot ^ 10*'"' erg/s. It is 
not clear whether these parameters can be achieved during the evolution [8, 9]. Finally, 
the numerical coefficient k is actually smaller than 1 [10]. For example, for a monopole 
magnetic field = 2'jTBar^{l — cos9) (and for Op = const) the exact expression (10) 
results in 



for a<^m, and 



(20) W,., = in/4-l/2)M}E_MBy,c 

for a = m (0.285 instead of 0.167). On the other hand, for a homogeneous magnetic field 
^ = irBy sin^ 9 one can obtain 



= 30 0| ImJ 

for a <C m, and 



(22) W^tot = (5/6 - 7r/4)M|_ziM52^2 ^ 

for a = m (0.0479 instead of 0.0333). Here we do not include into consideration the 
disturbance of the homogeneous magnetic field near a rotating black hole. This distur- 
bance is to diminish the magnetic flux passing through the hole i.e. to reduce the energy 
loss. Thus, in real objects a black hole may play a passive role only, and the main energy 
release can be connected with magnetic field lines passing through the accretion disk [11]. 
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Another criticism was connected with the vahdity of the Blandford-Znajek mechanism 
itself [12]. The point is that, as was demonstrated, during the derivation of expression 
(1), the condition (12) was actually used. But, clearly, the horizon is not in a casual 
connection with the external space and, hence, it cannot affect the flow structure or 
determine the energy flux flowing away from the rotating black hole. For the same 
reason, the surface currents cannot play any role in the black hole braking. As a result, 
the conclusion was drawn that a rotating black hole cannot work as a unipolar inductor. 

The goal of our paper is to clarify the ground of the Blandford-Znajek mechanism 
of the black hole energy loss. In particular, we arc going to determine the role of the 
horizon. In our opinion, a self-consistent analysis can be made on the ground of the Grad- 
Shafranov approach only. For this reason, in Sec. 2 we consider a simple example of the 
quasi-monopolc particle dominated accretion/ejection flow where the exact solution of 
the stream equation can be obtained. We will show that the natural boundary conditions 
at the pair creation region are enough to determine not only the longitudinal electric 
current /, but also the angular velocity of a flow flp as a solution to a problem. As 
a result, the energy loss of a rotating black hole is to be determined by the physical 
parameters at the particle creation region rather than the " boundary conditions" at the 
event horizon. Next, in Sec. 3 we consider some properties of a magnetically dominated 
flow when the flow structure is to be close to the force free one. Finally, in Sec. 4 we 
discuss the general properties of transonic flows. It is demonstrated that such flows can 
be realised only if there is no restriction of the longitudinal electric current in the source 
of plasma. 



2. — The particle dominated flow. Exact solution 

Let us consider the ideal magnetohydrodynamical cold flow in the vicinity of a slowly 
rotating black hole. For simplicity we shall consider the case when the energy density of 
the magnetic fleld is much higher than the plasma energy density 

(23) = %ii « 1. 

It does not mean that it is an electromagnetic energy flux that plays the main role in 

the black hole braking because for a nonrotating black hole the flux of electromagnetic 
energy vanishes. In this section we consider the particle dominated case when the energy 
flux of particles is much larger than the energy flux of the electromagnetic field: 

(24) ^2£^ » 1. 

^em 

The opposite case will be considered in Sec. 3. 

In the cold limit there arc four critical surfaces - two Alfvcnic and two fast magne- 
tosonic ones for ingoing and outgoing flows. Indeed, it is known that plasma can pass 
through the Alfvenic surface in one direction only - outward for an external Alfvenic 
surface and inward for an internal one (see e.g. [13]). Hence, plasma is to be created 
between two Alfvenic surfaces. As we shall see, it is the properties of the pair creation 
region that fully determine the flow structure including the energy loss of a rotating black 
hole. 

Unfortunately, the efficiency of pair creation in the magnetosphere of a black hole 
is not determined up to now. In particular, this process depends on the density and 
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energy of photons in a close vicinity of the black hole. As a result, if the density of hard 
gamma-quanta with energies £^ > IMeV is high enough, the particle creation can be 
connected with the direct process 7 + 7 + e~ [14]. On the other hand, if the density 
of hard photons is not high, the only possibility to create pairs is connected with the 
narrow sheet near the surface where the charge density 

(25) PC. = ^V. (^V^^- 

which is necessary to screen longitudinal electric field, changes the sign (for more details 
see [15, 16]). This surface is similar to the outer gap in the pulsar magnetosphere [17]. 
In what follows we consider the last mechanism of particle creation. 

The stream equation describing magnetic surfaces ^'(r, 0) in the vicinity of a rotating 
black hole was first formulated in [18]. In the 3 + 1 split language it can be written down 
in the compact form [19] 

(26) -V, I - ("f - ^)'^' - m2]v'^*1 + ^k^(Vvl/)2^ 
a [ avj^ J dw 

6471" 1 d fG\ 3 1 dry 



Here 



(27) A = a^-{^Y- ujfvj'^ - 

is an Alfvenic factor, 



(28) = 



2 47r77V 



n 



is an Alfvenic Mach number (n is the concentration in the comoving reference frame), 
and 

(29) G = a^zu^{E - n^Lf + a^M^L^ _ M^w'^{E - ujLf. 

In (26), the covariant operator acts in the three-dimensional metric (3), and the 
derivative d/d^ acts on the invariants i?(*), L{'^), and Of(4') only. 

Indeed, for a cold flow equation (26) contains four invariants. Two of them are the 
fluxes of the energy and angular momentum 

(30) E = E{'^) = ^+nr]{aj + u;vuu^), 

(31) L = i(*) = + nr]zuus, 

where n = {pm + P)/n is relativistic specific enthalpy. For the cold flow under consid- 
eration, n = const. The other two invariants are the angular velocity f2F(^) and the 
particle to magnetic flux ratio ??(4') 



(32) 



nup = r?(*)Bp. 
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As we see, in the general case the energy flux E{'^) contains not only the electromagnetic 
part Q.-pI/2'K, but the particle part iJ,r]{a'y + wzuu^) as well. The same takes place for the 
angular momentum On the other hand, according to (32), r] has different signs 

for ingoing and outgoing flows. 

Equation (26) is to be added by the Bernoulli one 

where 

(34) K = a^w^E - flpLf {A^ - M^) + [vj-'iE - wLf - a^L^] . 

This equation gives the implicit expression for the Mach number as a function of the 
magnetic flux * and four invariants 

(35) = M2[(V*)^ E{^), L(*), Of(*), ??(*)]• 

As a result, knowing the structure of the poloidal magnetic field and the four invariants, 
one can determine all the characteristics of a fiow. In particular [20] , 

^^g^ I _a^L-{nF-u)w^E-ojL) 



(37) 
(38) 



27r a2 - (fip - w)2w2 - M2 ' 
1 a^iE-flFL)- M'^{E-ujL) 

_ 1 {E - ^fL) (Of LM^ 



The general equation (26) is too complicated to be comprehensively analysed. Nev- 
ertheless, the solution can be obtained for a flow which is not too far from the known 
one. As a zero approximation one can take the (split) monopole magnetic field 

(39) * = *o(l-cos6l), 

which can be realised in the presence of the accreting disk terminating the ingoing and 
outgoing magnetic fluxes. It is clear that the monopole magnetic flcld (39) is an exact 
solution to equation (26) for a nonrotating black hole (and for E = const, i] = const, 
and Q.-P = L = Q). As was found earlier [21], for a monopole magnetic field and for 
rip = Oh/2 the pair creation region (i.e. the surface where the Goldreich- Julian charge 
density pcj (25) changes sign) is a sphere with a radius 

(40) nnj = 2VVh. 

On the other hand, using relations (33)-(34) for Of = Oh = 
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one can readily obtain the positions of Alfvenic and fast magnetosonic surfaces for an 
outgoing (M^ = 1) 



(42) 

and an ingoing (M^ = a^) 
(43) 



„(o>it) _ „(out) 



1/2 



87r2 ^E-^ - 



(in) 



a 



(in) 



1/2 



flows. Here the condition 7inj 1 and relation (23) resulting in 

(44) 
(45) 



(in) , 

ap <C 1 



were included into consideration. As we see, Alfvenic and fast magnetosonic surfaces for 
a nonrotating black hole have a spherical form and coinside each other. 
For slow rotation one can seek the solution of the full equation (26) as 



(46) 

Here 

(47) 



= <iQ[l- cose + fir, 6)]. 



m 



^0 is the total magnetic flux in the upper hemisphere, and / ~ 1. As a result, after the 
linearization of the general equation (26) we have 



(48) 



2 2 9 



J2 



d 



+ -2"o smd^— 



1 df 



de \sm.6de 



El 



+2 

— 2sign77 



[{2al - M^){flF - wa)^ - M^{flF - cof] sin^eicosei 



(fip — w)^ sin^ 9 cos 6 



(rip — uj){^Y — wa) sin^ ^cos^ 



-a; 



^5 



2 2 

ji r] 



rv-2„2„2 



2 ) sin^ ecose + al—M^ sin^ 9 cos 9 = 0. 



Here the values aoi^r), Moi^r), and Eq correspond to a nonrotating flow. Hence, in the 
limit e <C 1 it is possible to neglect the disturbance of the critical surfaces. Next, in 
(48) we have already used the expressions for the invariants L^™^ and which can 
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be obtained from the critical conditions on the Alfvenic surfaces. Indeed, analysing the 
nominators of the equation (36), one can obtain 

--)(out) (out) p 

(50) L(-) = - ^ „ *o sin2 e. 

It is interesting that within our approximation equation (48) does not depend on the 
disturbance of the energy SE = E — Eq. Finally, as wc shall sec, the critical conditions 
result in ftp const. For this reason in (48) we omit the terms containing dQp/dO. 

As equation (48) contains signry (to say nothing of the fact that the parameters 
cij^"*'' Ki and w^"'' s=; fJn arc different for ingoing and outgoing flows), wc actually 
have two different equations for internal and external regions. In other words, the whole 
equation (48) has a discontinuity at the pair creation region r = r-mj- Equation (48) 
contains all the information on the disturbance of the monopole magnetic field. 

In spite of the simplification, some properties of equation (48) are general and coinside 
with the properties of the full equation (26). At first, one can see that the region near 
the horizon m < r < / corresponds to the hyperbolic domain of equation (26). It 
means that, indeed, the horizon does not affect the magnetic field structure outside the 
black hole. On the other hand, it is clear that the number of boundary conditions does 
not depend on the simplification either. Using the relation [19] 

(51) 6 = 2 + i-s = 4 

for both ingoing and outgoing flows (6 is the number of boundary conditions, i = 4 is the 
number of invariants, and s = 2 is the number of critical surfaces), one can see that to 
determine fully all the characteristics of a flow it is necessary to specify eight boundary 
condition in the pair creation region. 

First of all, it is necessary to know the injection values of the particle concentrations 
and the Lorentz-factors of the flow. For simplicity, we shall further consider the case 



(52) n^- = ninj = const, 

(53) 7i^j = 7i„j = const. 



These four values determine the four invariants 

(54) ^(°"*) = ainj^ryinjTinj, 

(55) i;^'"^ = -ainjWinjTinj, 



(56) ?7^°"*^ = Vini = 

(57) »7^'"^ = -Vini 



Q:i„jni„j^7f„j - 1 

Q!inini„i^7?j - 



B 



H 



Here Bu = "^o/i^nr-^}. Bellow, we shall sometimes omit the symbol (out). As we can 
see, in the particle dominated case (24) the energies and have different signs 
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for ingoing and outgoing flows. It means that in this case it is the pair creation region 
that plays the role of the energy source. On the other hand, different signs of r]s represent 
the general property of the accretion/ejection magnctosphcrc. It is necessary to stress 
the difference from the model presented in [22] , where the velocity of plasma was assumed 
to vanish in the injection domain. 

Further, as a boundary condition one can use the absence of discontinuity of the 
magnetic flux 



(58) 



*(nnj-0) = *(ri„j+0) 



across the pair creation region r = rjnj . Then it is necessary to know two components of 
the surface electric current Jg flowing along the pair creation region: 



(59) 



Finally, it is necessary to know the potential drop along the magnetic field lines in 
this region 



(60) 



V{nni + 0) - y(ri„j - 0) = v„. 



The boundary conditions (52)-(53) and (58)-(60) are sufficient to determine all the 
properties of the accretion/ejection flow in the vicinity of a black hole. 

As has already been stressed, in the general case the solution of equation (26) cannot 
be obtained. On the other hand, the linearised equation (48) can be solved analytically. 
Indeed, separating the variables by the substitution 

(61) f{r,e)=g{r)sm^9cose, 

one can obtain the following ordinary differential equation for the radial function g{r): 



(62) 



2 2"^ 



w2N dff(r) 



a 2 fra 



+ 



an 



El 



4 {al - M2)2 El - altiW 



Ml 



El 



Ml 



) 9ir) 

(f^F -W)2 



^2 



1 



{al - MlY El alf,W 



(»F - io)^ 



a^jj^rf- 



"h 



2 2 



^0 



alii^r]'^ 



Ml = 0. 



The boundary conditions (58)-(60) for the injection radius r = rinj now have the form 

(63) 5(ri„j-0)=5(ri„j+0), 

(64) ^(n„j-0) = g(n„j+0) + Ai, 

(65) 7(ri„j - 0) = /(rinj + 0) + A/, 

(66) 0^,°"*) =0f)+A0F, 
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where the values Aj, A/, and Afip are to be determined from the mechanism of the pair 
creation. In particular, Afip is related to the potential drop in the pair creation region 

(67) A0f«0f-^" 



^nax 

Finally, these boundary conditions are to be supplemented by the regularity conditions 
on the fast magnetosonic surfaces Oq = Mq. 

In spite of the fact that the solution of equation (48) can now be easily obtained, wc 
are not going to discuss it in detail. The point is that the main properties of the flow 
for e <i; 1 do not actually depend on the radial function g{r). In other words, in the 
simple monopolc geometry (and for a slow rotation) the boundary conditions (65)-(66) 
together with the regularity conditions on the Alfvenic surfaces resulting in (49) and 
(50) are sufficient to determine the energy loss of a rotating black hole separately from 
the solution of equation (62). For this reason, we describe here the main properties of 
equation (62) only. 

First of all, one can see that equation (62) contains no singularity at the horizon. 
It is not siirprising because it results from the main property of the stream equation 
(26), namely, hyperbolisity near the horizon. Hence, it is not necessary to add any extra 
boundary condition for r = rn- On the other hand, it means that the disturbance of the 
monopole magnetic field remains small up to the very horizon 

(68) eV(rH,^)~£^<l. 

Secondly, in this approximation the positions of critical surfaces coinside with ones in zero 
approximation. Nevertheless, as was demonstrated, to obtain the solution it is necessary 
to use regularity conditions at Alfvenic and fast magnetosonic surfaces separately. As 
the Alfvenic singularity has already been used in (49) and (50), second, third, and fourth 
lines in (62) arc regular when Oq = Mq. Finally, at a large distance r^rx the solution 
of equation (48) docs not depend on the boundary conditions at all and coincides with 
the Bogovalov [23] solution in a flat space 



(69) e^f{r,9) 



2 — =- m — sm cos 0. 

\ c J 7ini VaJ 



A much more important thing is that our approach allows the determination of both 
the electric current / and the angular velocity fip and, hence, the energy loss of the 
rotating black hole. Indeed, using the expressions for angular momentum (49) and (50) 
together with the boundary conditions (65) (66), one can obtain 

/'7n^ o a;f^+a;i°"*^+AflF-si(2a;inj+A0F) ^ M<i - Ml^)AI 
(7U) S2f = — h 



2(1 -ei) al.M^-el)sm'0- 



Here 
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SO that £2 ~ £i. In the case under consideration, when the energy density of a secondary 
plasma is much smaller than the energy density of the poloidal magnetic field S2 <^ 1, 
we have w^"*^ <C Oh and w^"^ « Oh, so that 



(72) 
(73) 



Of = 



Oh + AOf + 



47r(a? 



2 

inj 



M2.)A/ 



r Op . 2 

4lT 



In particular, for A7 <C Iqj ~ OF^'o/47r and AOf Of we obtain 



(74) Of = Oh/2. 
Hence, according to (19), 

(75) ^^- = ^©'^>HC. 

Thus, one can conclude that the physical conditions at the pair creation region do allow 
the determination of not only the electric current flowing in the magnctosphere (this 
property holds in the flat space, sec e.g. [23]), but the angular velocity Of as well. 
Clearly, this conclusion is general and does not depend on our approximation. 

As we see, our results (73) (74) formally coinside with those obtained by Blandford 
and Znajek [1] within the force-free approximation. Nevertheless, they actually corre- 
spond to absolutely different tasks. Indeed, Blandford and Znajek found the values of 
current / = 2TrL and Of for which the magnetic field structure for a slowly rotating black 
hole in the force free approximation docs not differ strongly from the monopolc one. But 
within the force-free approximation, when we have only two critical surfaces (Alfvenic 
surfaces for ingoing and outgoing flows), in general case the current / and the angular 
velocity Op can be arbitrary. On the other hand, we have demonstrated that within the 
full MHD approach two extra critical (fast magnetosonic) surfaces fix the current / and 
angular velocity Op values for a double transonic flow. 



3. — The magnetically dominated flow 

In this section we consider some properties of the double transonic flows when the 
flux of the electromagnetic energy is much larger than the energy flux of plasma 

(76) ^ « 1. 

As is well-known, this relation can be rewritten in the form a ^ 1, where 

Oe'J'o 



(77) 



87rAmc^ 



is the Michel [24] magnetization parameter and A = n/nc] ^ 1. For simplicity, we 
again consider the slowly rotating black hole Oh^h <C 1. In this case the magnetospheric 
structure is to be close to the one obtained in the force-free approximation. Nevertheless, 
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there are two important differences. First of all, as previously, for any small but finite 
mass of particles the region in the vicinity of the horizon corresponds to the hyperbolic 
region of the general equation (26). On the other hand, the force free equation remains 
elliptical to the very horizon. Secondly, in the full MHD approach there are two additional 
critical surfaces - fast magnetosonic ones which are absent within the force-free approach. 

Unfortunately, no exact analytical solution can be obtained for magnetically dom- 
inated accretion/ejection in the vicinity of a rotating black hole. Nevertheless, it is 
possible to evaluate the main properties of the transonic flow by analysing the algebraic 
relations only. This approach has already been used in many papers (sec e.g. [25]). But 
in almost all of them the magnetic field structure has been considered as given. It does 
not allow a self-consistent analysis of the flow structure. As has recently been demon- 
strated in [26], some key properties, e.g. the position of a fast magnetosonic surface, 
cannot be determined correctly in a given (monopole) magnetic field. For this reason we 
now consider a more general case including the disturbance of the monopole magnetic 
field. 

So, let us consider again a flow with a monopole magnetic fleld in the pair creation 
region. Then we can rewrite the Bernoulli equation (33) in the form 



(78) - \ 



F 



g' 



1 2 /mV\ 1 f E-nFL Y_ 
+ 2" V i;2 j + 2 (fiF-a;)2tj7n E J ~ ' 

where we omit the term 5^ resulting in the unphysical root g <0. Here, by definition 

(79) g = 77^ 

(iZp — a;)^ti7^ 

and 

{E-uL)^ (Of-w)^w'(V*)2 



(80) e = 



GAn^E^ 

In particular, far from a black hole [26] we have 

2£^ df , , cos6' „ 
sm 6 06 sm 

so that ^ = for a monopole magnetic field. 

First of all, for an outgoing flow it is possible to use our results in a flat space [26]. 
Indeed, a fast magnetosonic surface corresponds to the intersection of two roots of equa- 
tion (78) at one point. On the other hand, equation (78) has two real positive roots if 
Q < 0, where Q is the discriminant of the third-order algebraic equation (78). Near an 
external fast magnetosonic surface r « rp, where the last term in (78) can be neglected, 
we have 

f82i n- ^ ^^"^^ ^ ^^"^^ (f+ ^ 



16 E'^ 432 E'^ V f2|r2sin' 
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the regularity conditions at the fast magnetosonic surface r = rp being 

(83) = 0, 

dQ/dr = 0, 
dQ/de = 0. 

As a result, we can rewrite the condition Q = near the external fast magnetosonic 
surface as 

(84) C(^F,g)+ ^,,A,2. =3^^^'^' 
and the condition dQ/dr = as 



02 r-2 sin^ 9 \E J 



(85) rpig) - , =0. 



dr 



As previously, the third condition dQ / 89 = is necessary for the determination of the 
magnetic disturbance. Using now the estimate r{d^/dr) ~ ^, one can obtain the position 
of the fast magnetosonic surface. For 7in <i; cr^/^ 

(86) rF{9) « i?LCT^/3 sin-^/3 9, 
when 9 > cr~^/'^, and 

(87) TF « i^L(a/7inj)'/^ 

near the axis. Here i?L = c/Op is the radius of the light cylinder. 

Further, as the root girp) for r = rp does not depend on the second term in (78), we 
have exactly 

(88) ,(,„,) 

Hence, 

(89) 7(^F, 0)=(—) ^'^ = T^'^ sin'/' 0, 

which corresponds to Michel's [24] result. The only difference is that this energy is 
achieved at a finite distance rp (86). As we see from (85), it takes place because we 
included the dependence of ^ on the field disturbance ef into consideration. Indeed, 
according to (84) and (85), it is the disturbance of the magnetic surfaces ef that plays 
the main role in (85) at the fast magnetosonic surface. Finally, according to (81), the 
disturbance e^f itself 



(90) eV(^F) « a 



-2/3 



is to be small at the fast magnetosonic surface. Moreover, as was found in [22, 26], 
outside the fast magnetosonic surface r rp, the particle acceleration and the magnetic 
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field coUimation become ineflFective. As a result, the disturbance of magnetic surfaces 
remains small up to infinity 



(91) 



-2/3 



< 1. 



Clearly, these results remain true for our problem as well. 

As to an ingoing flow, one can check that, according to (12), ^(rn,^) = 0. Hence, 
^ <C 1 for cr ^ 1 for an ingoing flow as well. As a result, using the same procedure as for 
an outgoing flow, where the discriminant Q now has a form 



(92p 



one can obtain for ^ 7^ 

(93) 
(94) 
(95) 
(96) 



2a' 



n 3 



432(f2F - a;)2ro2 E'' 

ff(^F) ~ -|, 

a{rF) 



(Of - w)2k72 



Here we introduce by definition e = E — VIyL. Together with (30)-(31) one can find that 

(97) i . ^. 

Thus, e/i? <C 1 for a magnetically dominated flow. 

First of all, we see that an internal fast magnetosonic surface for ^ 7^ is located 
much closer to the horizon than an Alfvenic one [25] 

(98) ^^{r^)=o?{rp:)'^, 

a 

so that o?{r-p) <C o?{rK). Here 

(99) o?(TK)^(Sl^-Q.vfvD\ 



corresponds to the position of the internal Alfvenic surface. Then, the Lorentz-factor at 
an internal fast magnetosonic surface differs only by the factor air-p) from its value in the 
injection region. But this fact has the coordinate reason only, and the formal increase of 
7 results from the difference in the ZAMO's positions in the regions of injection and fast 
magnetosonic surface. Thus, there is no intrinsic acceleration of ingoing particles, at any 
rate within the fast magnetosonic surface. Finally, according to (96), the disturbance of 
the monopole magnetic field remains small at the fast magnetosonic surface. Hence, as 
for the particle dominated fiow, the disturbance of the monopole magnetic field remains 
small up to the very horizon 



(100) 



£'/(r-H,^)«l. 
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But the most important property, in our opinion, is that the values of the invariants 
i(-), and 4°"*^ 

(101) « « 1^*0 sin^ e, 

(in) ^ o(°"*) ~ 



(102) fi^""^ « fip 



2 



(and, hence, the longitudinal electric current flowing in the magnetosphere) for a mag- 
netically dominated MHD flow are close to that considered by Blandford & Znajek for 
the force-free magnetosphere. It means that the double transonic flow is close to the 
force-free one. The difference occurs outside the fast magnetosonic surface only. But 
even in this region the disturbance of the monopole magnetic field remains small. 



4. — Discussion 

We have demonstrated that for a double transonic flow neither the longitudinal electric 
current /, nor the angular velocity is a free parameter, but is determined from the 
solution. Clearly, this result is general and does not depend on our approximation. On 
the other hand, exact values of the critical electric current I and the angular velocity 
flp corresponding to the double transonic flow certainly depend on the geometry of the 
magnetic field. Nevertheless, relations 



(103) 



fip ~ f^H/2, / ~ Igj, 



remain true for an arbitrary source of the external magnetic field. 

Unfortunately, for an arbitrary external magnetic field analytical calculations are 
impossible because, in particular, the position of the critical surfaces are unknown and 
they themselves are to be determined from the solution. Moreover, it is impossible to 
say that Alfvenic singularities determine the values of the angular momentum L and fast 
magnetosonic surfaces - the structure of the poloidal magnetic field. This was luckily 
realised for the monopole magnetic field only. In the general case, it is only all the critical 
surfaces taken together that determine invariants and the magnetic field structure. 

It is necessary to clarify the role of the "boundary condition" (12) which is formally 
used to determine the energy loss (1). This relation is true for any solution of the Grad- 
Shafranov equation which can be extended up to horizon [19]. The point is that relation 
(12) can be obtained by direct integration of the stream equation (26) which becomes 
parabolic at the horizon. Indeed, in the limit a — > we have 



(104) -Ve 
a 



a" 



i^(V*)2^ 
^ ^ d* 



327r^ d 



{E - nnLf 



0. 



As a result, multiplying (104) by 2^/(d^'/d^), one can obtain 



(105) 



2Aru'^{Q.T - Oh) 



dfip 



d4^ 
d^ 



+ 



d6i 



A^ 



= 0, 
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resulting in (12). 

Hence, in reality the "boundary condition" (12) c;ontains no additional information. 
The point is that wc do not know the magnetic flux 5'(rH.^) at the horizon (it is to 
be found as a solution of a problem), so that in general case relation (12) gives us no 
connection between the current I and the angular velocity fip. On the other hand, as was 
demonstrated, the necessary connection between / and fip results from the extra critical 
conditions on singular surfaces which are in casual contact with the pair creation region. 
The critical conditions give the same value of the longitudinal electric current because 
for a transonic inflow the disturbance of the monopole magnetic field at the horizon both 
for particle (68) and magnetically dominated (100) flows remains small. 

Finally, it is necessary to stress that our results depend sufficiently on the proposal 
that there is no additional restriction of the longitudinal electric current in the source of 
particles. In this case it is natural to assume the flow to be transonic and the current to 
be determined from the critical condition at the singular stirfaces. On the other hand, if 
the electric current / is determined by the particle creation process, the flow structure 
can be far from the transonic solution. 

This property is well-known in the flat space considered in connection with the pulsar 
magnetosphere. Indeed, if the electric current I is much smaller than the Goldreich- 
Julian one 

(106) 7gj = ^^osin^^, 

the force-free solution can be extended only up to the light surface |E| = |B| located in 
the vicinity of the light cylinder i?L = c/Q. [27]. A shock front must exist here resulting 
in a very effective particle acceleration. Simultaneously, the longitudinal electric currents 
/ flowing in the magnetosphere are closed in this region (in more details see [5]). As a 
result, outside the shock almost all the energy is transported by particles. 

On the other hand, if the electric current is larger than Jqj, there is a strong coUima- 
tion of the magnetic surfaces along the rotational axis [28] , the flow remaining subsonic. 
In this sense the Michel [29] force-free monopole solution with / = /gj terminates the 
flows with different asymptotics. Moreover, as has been demonstrated, a magnetically 
dominated MHD flow (cr 1) can pass through the fast magnetosonic surface only if 
the current I is close to the critical one. It is therefore not surprising that the transonic 
magnetically dominated inflow differs only slightly from the force-free Blandford-Znajek 
solution. In all other cases the properties of the magnetically dominated flow will be 
similar to the force-free one. 

Thus, one can conclude that within the MHD approach for a given current I ^ /qj 
there is an infinite number of solutions. Indeed, adding one extra boundary condition on 
the pair creation surface - the value of the electric current I -we simultaneously lose two 
critical conditions at the internal and external fast magnetosonic surfaces. In this case, 
the magnetically dominated flow will be similar to the force-free solution. In particular, 
for an ingoing flow it means that for small enough longitudinal electric currents the light 
surface |E| = |B| does not coinside with the horizon. Hence, in this case it is impossible 
to extend the solution up to horizon even within the MHD approach. 

In the table we classify different possibilities which can be realised if the longitudinal 
electric current is determined by the pair creation mechanism. For simplicity, we consider 
the case AJ = 0, AQp = 0, and the monopole geometry only. As we see, flows with 
small electric currents / do contain shocks in the vicinity of the light surface. On the 
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Table I. - Different possibilities of the flow with electric current fixed by the pair creation region 



KIgj 


/ = /gj 


/>/gj 




Of < 


Of < 


outflow with shock 


transonic outflow 


subsonic outflow 


inflow with shock 


inflow with shock 


arbitrary inflow 


Of = 


f^F = Qh/2 


fip = 


outflow with shock 


double transonic 


subsonic outflow 


inflow with shock 


flow 


subsonic inflow 


Qf > nH/2 


fiF > 


f^F > flH/2 


outflow with shock 


transonic outflow 


subsonic outflow 


arbitrary inflow 


subsonic inflow 


subsonic inflow 



other hand, flows with large enough currents are to be subsonic. Even for the current / 
coinciding with the critical one, the structure of the flow may be arbitrary. 

It is necessary to stress that there arc well known argiimcnts against stability of 
subsonic flows. Moreover, as has already been pointed out, the subsonic flow cannot 
pass the event horizon. Hence, either the subsonic inflow does not cross the horizon at 
all so that the magnetic field lines bend to the equatorial region [30], or this possibility 
corresponds to an unphysical solution with zero velocity and infinite particle density at 
the horizon [13]. It is clear that only the solution of the stream equation (26) could 
determine the real structure of the magnetic field in this case. Unfortunately, up to now 
such solutions have not been constructed. On the other hand, solutions with shock (i.e., 
those containing light surfaces |E| = |B|) are real and can be realised. Then, outside the 
light surface the Grad Shafranov approach itself becomes invalid. These solutions arc 
interesting in connection with the possibility of an effective particle acceleration near the 
light surface (see e.g. [30, 31]). 

Incidentally, the fact that there is a restriction of the magnetic field structure in the 
vicinity of the horizon has already been discussed by Hirotani at al [25] for a magnetically 
dominated flow. In their interpretation, the magnetic field cannot have an arbitrary 
structure (in particular, homogeneous) near the horizon of a black hole. In our opinion, 
their results actually confirm our point of view - in the presence of the light surface 
|E| = |B| it is impossible to prolong the MHD solution to the very horizon. But it does 
not mean that there is a restriction of the magnetic field structure near the horizon. 

As to the restriction found in [25], it can be obtained by an even simpler way. As 
is well known, on the horizon the electric field is equal to the magnetic one for any 
geometry of the magnetic field: it results from the membrane paradigm. On the other 
hand, one can extend the solution up to horizon only if the electric field is lower than 
the magnetic one. Hence, one can write down this condition as 

(107) ^{B^ - E^)r^ > 0. 

Using now the definitions (6) and (7), one can rewrite (107) in the form 

- 4^ sin^cos^P + 2^ - 2 f 1 - 



(108) 
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I /"^^'^ 2 ( -2^ ^ +2^>0 



where the parameter 

(109) P = -rn 



depends on the magnetic field structure near horizon. One can easily check that the 
condition (108) coincides identically with the one obtained by Hirotani et al [25]. 



5. — Conclusions 



Thus, it was demonstrated that: 

1. It is impossible to consider relation (12) as a boimdary condition. This relation is 
true for any solution of the Grad-Shafranov equation which can be extended up to 
horizon. In reality, we do not know the magnetic flux ^'(rH,^) at the horizon (it 

is to be found as a solution of a problem), so that in the general case relation (12) 
gives us no connection between the current / and the angular velocity flp. 

2. For the finite mass of particles in the very vicinity of the horizon there is a hy- 
perbolic region of the stream equation which is absent altogether in the force-free 
approximation. Hence, the stream equation needs no boundary condition at the 
horizon. 

3. It is impossible to say that it is the surface electric current that results in the 
braking of the rotating black hole. The surface electric currents Jh (14) as well as 
the surface charges cth (13) are unphysical and cannot act on the surface of a black 
hole. They only give a convenient way to describe the flux of the negative energy 

falling onto the black hole horizon. 

Hence, the Blandford-Znajek mechanism of the electromagnetic energy extraction 
from the rotating black hole faces no problem connecting with the causality disconnection 
between the event horizon and the outer magnctosphcrc. In particular, for a double 
transonic flow the values of the longitudinal electric current / and the angular velocity 
Clp (and, hence, the energy loss Wtot) are to be determined by the physical parameters 
in the pair creation region and by the critical conditions at the singular surfaces which 
are in casual contact to each other. 

On the other hand, if the longitudinal electric current is determined by the pair 
creation region, then there is an infinite set of solutions which do not pass one or both 
fast magnetosonic surfaces. In particular, for small enough longitudinal electric currents 
the effective particle acceleration can take place in the vicinity of the light surfaces. In 
this case the ideal MHD solution cannot be extended to the horizon. 

As in the pulsar magnetospherc, if the secondary plasma is enough to screen the 
longitudinal electric field, its charge density and electric currents produce the fiux of 
the electromagnetic energy propagating from the central star to infinity. For the same 
reason, a rotating black hole embedded into an external magnetic field works as a unipolar 
inductor extracting its energy of rotation by the flux of the electromagnetic energy. 
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